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und diese beiden Gleichungen lassen sich zu einer vereinigcn, 
wenn wir die erste mit * multipliciren und zur zweiten 
addiren: 

(2) c curl (g + m) = i 8 ( <£ + * m > , 

Figure 1. The first appearance in print of the Riemann-Silberstein vector 



1. Introduction 

This review is devoted to various applications in classical and in quantum theory of 
the following complex combination of the electromagnetic field vectors: 

F = ^r + [ ^ ^ 

where e is the dielectric constant and (i is the magnetic permeability. The 
denominators in the definition of the RS vector are needed to make the dimensions of 
both terms equal. The RS vector is most useful when dealing with pure radiation field 
in empty space or in a homogeneous and static medium where e and fi are constant. 
In this review we restrict ourselves most of the time to this case. 

The RS vector has been studied in some detail for the first time in 1907 by 
Ludwik Silberstein [1, 2]. In the second paper Silberstein wrote that this vector 
made its first appearance (see Fig. 1) in the lectures on partial differential equations 
by Bernhard Ricmann edited and published in 1901 by Heinrich Weber [3]. Later 
Silberstein incorporated this vector into the quaternionic formulation of the relativity 
theory [4]. The name the Riemann-Silberstein (RS) vector was introduced by us in 
[5] and this name is now generally accepted. 

This review is divided into two parts: classical theory of electromagnctism and 
quantum electrodynamics. Many results presented here are not original — most of 
them are taken from our earlier publications [5-23] . We shall show that the RS vector 
offers a uniform description linking the wave aspect and the corpuscular aspect of 
electromagnetism. This link has also been emphasized in a recent article by Smith 
and Raymer [24]. We must also mention here very extensive reviews by Keller [25] 
and Saari [26] where the RS vector is used to describe quantum mechanics of photons 
with special emphasis on photon localization. 

The Maxwell equations in the absence of sources read: 

d t D{r,t) =VxH(r,t), (1.2a) 

d t B(r,t) =-VxB(r,i), (1.26) 

V-D(r,t)=0, (1.2c) 

V-B(r,t)=0. (1.2d) 

By adding the first equation multiplied by l/\/2e to the second equation multiplied by 
i/-\/2fi and doing the same for the third and fourth equation one obtains two complex 



equations 

id t F(r,t) =VxG(r,t), (1.3a) 

V-F(r,i)=0, (1.36) 
where the second complex vector is defined as: 

G=^=+i^=. (1.4) 

In vacuum (or in a homogeneous and static medium), due the simple form of the 
constitutive relations, 

H = B/ft, E = D/e, (1.5) 

the vector G is a multiple of F (G = cF). In this case the four Maxwell equations 
are equivalent to the following equations for the complex RS vector: 

id t F(r,t) = cV x F(r,t), (1.6) 
V-F(r,t)=0. (1.7) 

Thus, with the use of the RS vector, one reduces the number of equations by a factor 
of two at a price of working with a complex field. Note that when e and \i vary in 
space or in time the reduction to two equations does not work [5] . 

The use of a complex vector to describe the electromagnetic field is not a 
complication since it is often customary, as a matter of convenience, to seek complex 
solutions of Maxwell equations. The imaginary part is rejected at the end as 
nonphysical. We are not so wasteful when we are working with the RS vector since both 
parts are physically meaningful. The real part of the RS vector describes the electric 
field and the imaginary part describes the magnetic field. The economy, however, is 
not the only reason to use the RS vector. 

The usefulness of the RS vector stems from several factors. 

• The use of such a complex representation of the electromagnetic field greatly 
simplifies the derivation of various solutions of Maxwell equations. This has been 
exploited already long time ago by Bateman [27]. 

• The complexification of the electromagnetic field allows for generation of new 
solutions from the old ones. 

• The use of the RS vector underscores the role of the Hertz and Debye 
superpotentials. 

• The RS vector is as close to the photon wave function in coordinate space as 
possible, taking into account the limitations imposed by the lack of the position 
operator for photons. 

• The quantization of the electromagnetic field is more transparent when the field 
operators are combined into the RS field operator. 

In this review we will illustrate the merits of the RS vector with many specific 
examples. 



2. Basic properties of the RS vector 



In the absence of sources the splitting of the electromagnetic field into the electric and 
magnetic part is not absolute because in this case Maxwell equations are invariant 



under the duality rotation of the electromagnetic field ([28], [29] p. 274), 

E' = E cos <p — cB sin <p, (2.1a) 

cB' = cB cosy + E sin (p. (2- Id) 

The duality rotation expressed in terms of the RS vector means the multiplication of 
this vector by a constant phase factor, 

F -> exp(itp)F. (2.2) 

The duality rotation changes the relative proportions of the electric and magnetic 
field. The presence of electric charges in the real world breaks the symmetry under 
the duality rotations and introduces the physical distinction between electricity and 
magnetism. In other words, in the duality plane where the dual rotations (2.1a) take 
place there is a distinguished direction introduced by the electric charge. 

There is some analogy between the change of phase of the RS vector and the 
change of phase of the wave function in wave mechanics. Namely, the most important 
physical quantities: the field energy (Hamiltonian) H, the momentum P, and the 
angular momentum M are invariant under the change (2.2) of the phase of F. Of 
course, the phase of F, unlike the phase of the wave function, is observable once we 
distinguish between the electric and magnetic field. 

The quantities H, _P, and M are the generators of translations in time and 
space and rotations. The full list of generators includes also the generator of Lorentz 
transformation TV. All the generators are built from the energy density, 

D 2 B 2 

and the momentum density (Poynting vector divided by c 2 ), 

(D x B) = ~(F* x F), (2.4) 

1C 

as the following space integrals: 

H = Jd 3 rF*-F, (2.5a) 

P=- [d 3 r(F* xF), (2.5b) 
ic J 

M=- d 3 rrx(F*xF), (2.5c) 
ic J 

N= d 3 rrF*-F. (2.5d) 



The generator of Lorentz transformation TV has the simple form of the first moment 
of energy only at t = 0. Since it is not a constant of motion, at other times it develops 
a time-dependent part —c 2 tP. 

The square of the RS vector is a complex combination of the two invariants S 
and P of the electromagnetic field: 

F 2 = S + iP = e [\ {E 2 - c 2 B 2 ) + mE ■ B] . (2.6) 

Therefore, it is not surprising that the transformation properties of the electromagnetic 
field expressed in terms of the RS vector have the form of rotations preserving the 
square of F . Under ordinary rotations the RS vector transforms as a normal vector: 

F' = F cos tp + n x Fsinifi + n (n ■ F) (1 — cosy), (2.7) 



where tp is the rotation angle and n is the unit vector of the rotation axis. A pure 
Lorentz transformation is obtained from this formula by substituting for the rotation 
angle tp the angle of an "imaginary" rotation —hp: 

F' = Fcosh^ - in x F sinh tp + n(n ■ F) (1 - coshV), (2.8) 

where ip = arctanh(v/c) and n is the unit vector in the direction of the velocity. 



Classical theory 



3. Solution of Maxwell equations by Fourier transformation 



The RS vector is ideally suited to use the method of Fourier transformations. Since the 
RS vector is complex- valued one does not have to impose extra reality conditions on its 
Fourier transform F(k, t), as in the case with electric and magnetic fields. To simplify 
notation we use the same symbol, but a different argument, to denote the Fourier 
transform. Taking the three-dimensional Fourier transform of both sides of Maxwell 
equations one arrives at a simple ordinary differential equation and the transversality 
condition: 



(3.1) 



d t F(k,t) =ckxF(k,t), 
k ■ F(k,t) = 0. 

The transversality condition requires that the vector F(k,t) lies in the subspace 
orthogonal to k. We may build the basis in this subspace from two orthogonal real 
polarization vectors, but it will be more convenient to choose a complex basis e±(fe) 
in this subspace built from the eigenvectors of the cross-product, 

k x e±(fe) = =pifce±(fc), (3.2) 

All normalized vectors that satisfy these equations differ only by 
We choose the phases of these vectors such that e+(fe) = e(fc) and 
= e(— k). The e(k) vector will be chosen in the form: 

k x (m x k) — ik(m x fe) 

e(«) = TkT, r, . ( 3 - 3 ) 

V2 k\m x k\ 

Choosing the m vector along the z axis, one 



where k = \k 
phase factors. 
e_(fc) = e*(fe) 



where m is an arbitrary unit vector 
obtains: 
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(3.4) 



The polarization vector e(fe) is dimcnsionless and is normalized to one, e*(fe)-e(fe) = 1. 



In spherical coordinates k x 
zation vector has the form: 



k cosy sin 9, k y — ksinipsmO, k z = kcos0, the polari- 



e(M 
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COS ip COS 6 

sin tp cos 9 
sin ( 



- i sin tp 
icos tp 



(3.5) 



Of course, this choice of e(fc) is not unique because it can always be multiplied by an 
arbitrary fc-dependent phase factor. 



After decomposing F(k, t) in the basis of polarization vectors, we solve the 
resulting ordinary differential equations and we arrive at: 

F(k,t) = e(k)f + (k)e- iut + e(-k)f*_(-k)e iut , (3.6) 

where u = ck. The use of complex conjugation and the sign reversal in the second 
term may seem awkward, but this notation will be justified later in quantum theory. 
The inversion of the Fourier transformation gives the following integral formula of the 
RS vector representing the general solution of the Maxwell equations: 

/H 3 k 
j^y^e(k) [/+(fe)e lfe — « + /* (fc) e - Ifc -+*"] . (3.7) 

In the second term we changed the integration variables from k to —k. This 
integral representation is an expansion of the RS vector into circularly polarized 
monochromatic plane waves, 

F xk (r,t) = e(k)e ix ^ r - ut \ (3.8) 

where A = ±1 specifies two circular polarization states. The appearance of only one 
polarization vector (and not two) in the Fourier representation (3.7) is a characteristic 
feature of the RS vector [22] not shared by the real electromagnetic field vectors D 
and B. This property simplifies many calculations. 

Two complex amplitudes f± represent the independent degrees of freedom of the 
electromagnetic field. They contain all information about the field. In particular, 
the dynamical quantities (2.5a-d) have the following representation in terms of the 
functions f\(k): 



ff = E/T /;;(fe)u/A(fe) ' (3 - 9a) 

P = E —mk)kfx(k) 7 (3.96) 



d 3 k 



A 



/rl''k 
— f*(k) (iD x xk + Xn) f x (k), (3.9c) 

/r] 3 k 
Y-fl(k)uiD x f x (k), (3.9d) 



A 

where n is the unit vector in the k direction. The symbol D x denotes the covariant 
derivative on the light cone [6, 19], 

D x = d i\a(k), (3.10) 
a (k) = ^emOe.m = ^^, (3.1D 

i 

where d denotes the nabla in the k space. 

To connect the representation (3.8) of a plane wave expressed in terms of the 
RS vector with the formulas found in standard textbooks [29] we choose the direction 
of propagation along the z axis which has already been distinguished in the formula 
(3.4) for the polarization vector. The simplest choice of e(fc) in this case gives the RS 
vector in the form: 

" 1 



Separation of this expression into its real and imaginary parts leads to the standard 
formulas for the electric and magnetic fields of the circularly polarized monochromatic 
plane waves (we dropped the overall factor), 



E(r,t) 



COs(fcz — Lot) 

-A sin(fcz — Lot) 




B(r,t) 



A sin(fcz 
cos(fcz - 




-LOt) 
LOt) 



(3.13) 



One obtains the left-handed or the right-handed circularly polarized waves propagating 
in the z direction depending on whether A is positive or negative. 

4. Superpotentials 

It has been known for more than one hundred years that the use of potentials simplifies 
the process of solving the Maxwell equations. A systematic method to find and classify 
all solutions of the Maxwell equations with the use of the RS vector can be based on 
just one complex scalar function. 



4-1- Hertz potentials 

We introduce first a complex vector supcrpotcntial (complex Hertz potential) II. The 
real and imaginary parts of II can be identified with the electric and magnetic Hertz 
potentials II e and II m [30, 31]. The RS vector is constructed from II according to 
the formula: 



F(r,t) 



Vx 



(VxII(r,t)). 



(4.1) 



The Maxwell equations obeyed by F impose a condition on II. In order to obtain this 
condition we act with the differential operator [l/cd t + iVx] on both sides of (4.1). 
The left hand side vanishes on account of (1.6) and we obtain the d'Alembert wave 
equation for the curl of II, 

(Vxn(r,f))=0. (4.2) 



A 



The Hertz potentials deserve the name "superpotentials" because the electric and 
magnetic field vectors arc their second derivatives and not the first derivatives as in 
the case of standard potentials. 

There is a lot of freedom in choosing the complex vector II. The class of gauge 
transformations of II that leave the RS vector unchanged is very large. We shall 
reduce this freedom by requiring that the vector II is derived from a single complex 
function satisfying the d'Alembert equation. The simplest choice is to assume that II 
has some fixed direction, II = m<I>. Then, 



F(r,t) 



-d t (V x m) + V x (V x m) 



*(r,t). 



(4.3) 



Choosing the direction of m along the z axis, as in the formula (3.4), one obtains: 

d x d z + i/edtdy 

-i/cd t d x $(r,t). (4.4) 
-dl - dl 

For every complex function $ satisfying the scalar d'Alembert equation, 



F(r,t) 



d y d z 



AU(r,i) = 0, 



(4.5) 



the RS vector (4.4) satisfies the Maxwell equations (1.6). The formula (4.4) is a 
complexified version of the representation introduced a long time ago by Whittaker 
[32, 33] so that the real and imaginary parts of $ could be called the electric and 
magnetic Whittaker potentials. The Whittaker construction is particularly useful in 
the description of beams of electromagnetic radiation since then there is a preferred 
direction of propagation. The potential $ satisfying the d'Alembert equation can be 
expanded into plane waves, 



Differentiating both sides of this formula according to (4.4) we reconstruct the Fourier 
representation (3.7) of the RS vector. Therefore, every solution can be obtained by 
the Whittaker construction from an appropriate scalar solution of the d'Alembert 
equation. The two amplitudes f±(k) in the Fourier expansion of $ are the amplitudes 
of circularly polarized parts of the electromagnetic field. 

4-2. Debye potentials 

There is also another choice which is very convenient when dealing with radiation 
produced by localized sources. In this case one assumes that the superpotential II has 
the radial direction: II = r<£>. It follows from (4.2) that $ must satisfy the d'Alembert 
equation (4.5), as in the Whittaker construction. The real and imaginary parts of $ 
are known as the electric and magnetic Debye potentials [34]. The RS vector in this 
case has the form: 



The advantages of using scalar functions to represent the electromagnetic field 
are obvious. Scalar functions represent the independent degrees of freedom. Moreover, 
solving the d'Alembert equations for a single function is much simpler than solving 
several coupled equations. The use of the RS vector results in the replacement 
of two real scalar functions, electric and magnetic superpotentials, by one complex 
combination. 

5. Assorted solutions in the new framework 

We shall illustrate now the advantage of using the RS vector by constructing 
various important solutions of the Maxwell equations with the use of superpotentials 
constructed from the scalar solutions of the d'Alembert equation. We will be rather 
careless about the normalization of the RS vector since, due to the linearity of Maxwell 
equations, the overall complex amplitude is arbitrary. Obviously, from a given solution 
i> of the d'Alembert equation one may construct two different solutions of Maxwell 
equations by using either the Whittaker or the Debye construction. 

A practical method of finding explicit solutions of the d'Alembert equation is the 
reduction to ordinary differential equations by the separation of variables. There are 
368 coordinate systems in flat space [35] that allow for the separation of variables 
in the d'Alembert equation. From this vast number of solutions we will select a few 
examples. 




(4.6) 



F(r,t)= -dt-iVx L$(r,t), 
l c 

where L is the angular momentum operator (up to a factor H/i) 
L = r x V. 



(4.7) 



(4.8) 



5.1. Monochromatic waves 



Monochromatic waves are obtained by the separation of the time variable in the 
d'Alcmbcrt equation. Thus, we seek the solution in the form: 

$(r,i) = *(r)T(t). (5.1) 

The substitution of this product into the d'Alembert equation results in the following 
two separate equations: 

d 2 t T{t) = -u; 2 T(t), (5.2a) 
(A + lu 2 /c 2 ) tf(r) = 0, (5.26) 

where lo 2 is a separation constant. To avoid an unphysical exponential growth of 
T(t), we assume that u is real. The first equation for a given uj 2 has two independent 
solutions: exp(=Fiwi), hence $(r,t) has two terms: 

$(r, *) = *+(r)c- itJt + *_(r)c i " t . (5.3) 

The RS vector obtained from such a superpotential (by either Whittaker or Debyc 
construction) has the form: 

F(r,t) = F+(r)e- iu]t + F_(r)e iu '. (5.4) 

It follows from Maxwell equations that F± are eigenfunctions of the curl operator, 
known also as the Trkalian fields [36] , 

cV x F±(r) = ±u)F±(r). (5.5) 

The real part and the imaginary part are both Trkalian fields. 

Monochromatic solutions of Maxwell equations have an interesting property that 
follows immediately from (5.4). The electric and magnetic field vectors at each space 
point trace ellipses defined by the following parametric representations: 

Electric : 5i(F + + F_) cos(o;i) + 3(F+ - F_) sin(wi), (5.6a) 

Magnetic: $X(F+ + F-) cos(cjt) - - F_) sin(wt). (5.66) 

In general, electric and magnetic ellipses are different, but they coincide when there 
is only one term present in (5.4). 

One can further refine the classification of monochromatic waves. One may call a 
monochromatic solution of Maxwell equations stationary if its RS vector has the form 
F(r) cxp(— iXut). Each monochromatic wave is a linear combination of two stationary 
waves as shown in (5.4). Stationary waves are direct counterparts of stationary states 
(energy eigenstates) in quantum mechanics. The energy density and the Poynting 
vector are time- independent for a stationary wave while in general for a monochromatic 
wave they oscillate with twice the wave frequency. 

The Hclmholtz equation (5.26) for ty(r) can be separated in eleven coordinate 
systems [37]. Of course, the simplest case is the separation in Cartesian coordinates 
that leads to the solutions of the d'Alembert equation in the form of monochromatic 
plane waves (3.8). 

5.1.1. Bessel beams. Electromagnetic waves described by Bessel functions were 
known in the context of electromagnetic waves inside the cylindrical wave guides 
[38]. Their use in the description of light beams came later [39]. Bessel beams are 
monochromatic and non-diffractivc. The mathematical description of Bessel beams 



is obtained by separation of variables of the Hclmholtz equation (5.26) in cylindrical 
coordinates. This leads to the following solution of the d'Alembert equation: 



z—<jjt-\-Mip) 



(5.7) 



k 2 and 



where fey is the component of k in the direction of propagation, k±_ = yk 2 

Jm(k±p) is the Besscl function. We have chosen the normalization factor in such a 
way that in the limit, when k± — > one obtains monochromatic plane waves (3.13). 
The parameters ku,k±, and M are the analogs of quantum numbers. We shall give a 
precise formulation of this analogy later when we describe photon states. The meaning 
of fen, kj_, M and A will then become clear. 

The Whittakcr construction leads to the following expression for the RS vector 
of the Bessel beam: 



F, 

Fy 



^iA(A;||2— iut-\-Mtp) 



K-e ltp J M +\(kj_p) + n + c- lv J M - X (k±p) 
-in„e lip J M+ x{k ± p) + iK + e" lip JM-\{k±p) 
-in ± /kJ M (k±p) 



(5.8) 



where k± = (1 ± fcii/fc)/2. These formulas slightly simplify when one uses the 
projections (F p , F V ,F Z ) of F on the basis vectors of cylindrical coordinates (e r , e v ,e z ), 



F n 



F z 



■ 3 iA(fc|i z— ut-\-M (p) 



(5.9) 



K-J M +\(kxp) + K + JM-\(kxP) 

-in^J M +\(kj_p) + iK + J M -\(k±p) 
-ik±/kJ M (kj_p) 

In the limit k± — > 0, only the Bessel function Jo survives. Choosing M = 1, A = 1 
or M = -1, A = -1 one obtains circularly polarized plane waves (3.12) propagating 
in the z direction. Bessel beams cannot be exactly realized in Nature since, like plain 
waves, they carry infinite energy flux — Bessel functions fall off too slowly when p — > oo 
to make the integral of the Poynting vector over the plane perpendicular to the beam 
axis finite. However, beams of light resembling Bessel beams in the neighborhood of 
the propagation axis are now commonly produced in optical laboratories [40] . 



5.1.2. Multipole fields. Multipole fields are very useful in the studies of electro- 
magnetic radiation by localized sources. The derivation of the solutions of Maxwell 
equations in the form of multipole fields is usually a cumbersome task (cf., for example, 
[29]). The use of the RS vector expressed in terms of the complex Debye potential 
significantly simplifies the calculations. All one has to do is to choose the solution of 
the d'Alembert equation $ in the form: 

- iXut Y lm (6,v)ji(kr), 



(5.10) 



where Yi m is the spherical harmonic and ji(kr) is the spherical Bessel function and 
then to perform the differentiations required in evaluating F. Treating the expression 
(5.10) as the complex Debye potential in the formula (4.7) we obtain: 



F(r,t) 



iXujt 



(AHVx; 



L., 
L- 



Y lm (e,ip)ji{kr) 



(5.11) 



These are the multipole solutions of Maxwell equations. In the description in terms 
of the RS vector the existence of two types of multipoles: the electric type and the 
magnetic type finds its counterpart in the two values of A. The formula (5.11) has 
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F e 




r 



&(r,0,<P,t)- 



(5.12) 



a simple form when expressed in terms of the projections (F r ,F0,F v ) of F on the 
spherical coordinate versors (e r , eg, e v ), 

-1(1 + 1) 
—Xkr/s'mOdp — dgd r r 
Xkrdg — 1/ sin8 d v d r r 

The function (5.10) is well known from wave mechanics where it represents the 
solution of the time-independent Schrodinger equation in free space with quantum 
numbers fc, I, and m. We have chosen the standing waves described by the spherical 
Bessel functions, but one may also replace them by spherical Hankel functions 
(kr) to describe the outgoing or incoming radiation. 



5.2. Nonmonochromatic beams 

Monochromatic solutions of Maxwell equations are a very convenient theoretical tool, 
but they are unrealistic since they describe a never ending wave train. In many cases 
the idealization involved here is harmless, but quite often nonmonochromatic solutions 
of Maxwell equations are more appropriate. The time-dependence of such solutions is 
nontrivial and the RS vector comes in handy. 



5.2.1. General plane waves. The simplest example of a nonmonochromatic wave is a 
general plane wave. Such a wave F(k-r / c—t) propagating in a given direction k = k/k 
is obtained from the general solution (3.7) by leaving out the angular integration and 
retaining only the integration over uj = fee, 

/>oo 

F(k-r/c-t) = e(k) duo A + (uj)e iuj ^ r/c -^ + A-{(j)e-' m ^ r/c -^ , (5.13) 



where, as usual, the two terms correspond to two circular polarizations. We placed 
the polarization vector e(k) outside the integral because it does not depend on k. By 
allowing lo to take on positive and negative values, we may write the result as a single 
integral, 

/oo 
duA^y^'^ = e(k)A(k ■ r/c - t), (5.14) 
-oo 

where the wave amplitude A is the Fourier transform of the spectral function A(u). 



5.2.2. Gaussian beams. Gaussian beams, unlike Bessel beams, are neither mono- 
chromatic nor non-diffractive, but they are more realistic; they carry finite energy 
flux. Gaussian beams are used to describe light emitted by lasers. Very often they 
are treated as monochromatic, but this is only possible in the paraxial approximation. 
Here we consider Gaussian beams as exact solutions of Maxwell equations. The RS 
vector is ideally suited to describe these solutions. 

The solutions of Maxwell equations that are not monochromatic form a more 
diverse collection because one is then not restricted to a rather rigid framework 
imposed by the Helmholtz equation (5.2 b). Gaussian beams can be obtained after 
the change of variables, t = z + ct, £ = z — ct, in the d'Alembert equation. This leads 
to: 



(4d T d ( + A x )$(r ± ,T,0 = 0- 
where r±_ = (x, y) and Aj_ is the Laplacian in two-dimensions. 



(5.15) 



We will seek the solutions of (5.15) by separating the dependence on £, 

10^(0 = - M 0(C), (5.16a) 

id T ip(rj_,T) = -— A_L^(r_L,r). (5.166) 

We have chosen a real value of the separation constant p (it can be positive or negative) 
because any imaginary addition will cause the unphysical exponential growth of the 
field. The resulting Schrodingcr-likc equation in two dimensions (5.166) can be solved 
by the separation of variables in 26 coordinate systems [41]. Gaussian solutions, 
however, can be obtained without resorting to variable separation. To this end we 
express the solution of (5.166) as a Fourier transform and obtain the following formula 
for $: 

$(r±,z,t) = e^ {z - ct) Jd 2 kf{k ± )e ik ^ r ^e X p(^-i^{z + ct)^, (5.17) 

where f(k±) is an arbitrary function. To obtain a simple Gaussian beam we choose 
/ as a pure exponential without polynomial prefactors (such polynomials can be 
generated later by taking derivatives with respect to r±), 

f(k x ) = e- k± -x- k± . (5.18) 

The complex 2x2 matrix K is arbitrary but its real part should be positive to produce 
a convergent integral. With this choice of f(kj_) the Fourier integral can be evaluated 
and the function $(rj_, z, t) takes on the form (dropping an overall numerical factor): 

gi^(z-ct) / 1 \ 

*(r Xj z, t) = cxp -r ± • — - r ± \. (5.19) 

^Det {aK + Kz + ct)/pi} V 4K + i{z + ct)/vI ) 

The simplest Gaussian beam is obtained when K is a multiple of the unit matrix, 
K = (l 2 /4)I. Then, 

^^*) = ^^ eXP (-^)j' (5 - 20) 

where h(z + ct) = I 2 + \(z + ct)/fi. The RS vector constructed from this function is: 

-ipe-^^/i 2 (r) + ipe^(2/i(r) - p 2 )/fi 
pc-^fih 2 ^) + pc ilfi (2h{T) - p 2 )/fi 
2h( T )(h(r) - p 2 ) 

This compact form shows the efficiency of the description in terms of the RS vector. 
The formulas for the electric and magnetic field would be much more complicated. 
However, in order to calculate the energy density and the Poynting vector one does 
not need the real and imaginary parts; one may use directly the formulas (2.3) and 
(2.4). We will present here the results for the energy per unit length along the z axis 
and the energy flux through the xy plane. Keeping the normalization the same as in 
(5.21) one obtains: 

dxd yF<.F = , 2l ^ + 1 f/ + 3 , !U f t xF = ff ^-3, (5 . 22) 



exp(i/iC- 7TFT 

F(^,C,r)=2 \ 4 . ( ' 

h 4 (r) 



(5.21) 



Pp? ' i J a l 8 p 2 

These values do not depend on z or t as is required by the continuity equation for 
the energy and the energy flux. Standard Gaussian beams obtained in the paraxial 
approximation do not have this property. Note that the beam does not transmit any 
energy when Ip — (3/2) 1 / 4 . 



5.3. Construction of electromagnetic beams as superpositions of simpler solutions 

Monochromatic plane waves, Bessel beams, Gaussian beams, multipolc fields arc 
examples of complete sets of solutions: every solution of Maxwell equations can be 
constructed as a superposition of solutions from a given set. Construction of solutions 
in this way is most easily done with the use of the RS vector because the building 
blocks are most often complex functions and postponing the separation of the electric 
and magnetic parts until the very end greatly simplifies the calculations. 

Monochromatic plane waves are the most natural building blocks owing to their 
very simple mathematical form and their close relation to the Fourier transform. Every 
solution of Maxwell equations can be built from plane waves according to the formula 
(3.7). An expansion into monochromatic plane waves is the best source of information 
about the beam structure that helps in their construction. It also gives the information 
about the beam spectrum. Usually one cannot evaluate the integrals over k in a closed 
form. We shall present here four examples where the integrals over k can be explicitly 
evaluated. Examples of solutions obtained as superpositions of solutions other than 
plane waves will also be given. 

5.3.1. Construction of Bessel beams from plane waves. Bessel beams are monochro- 
matic and the component of the wave vector fcii along the beam axis is fixed. Therefore, 
there is no integration over k or kn . This means that k± is also fixed. The only 
remaining integration variable in (3.7) is the angle <fi. The simplest way to construct 
Bessel beams from plane waves is to start from the following integral representation 
of the Bessel function: 



This formula can be extended to the three-dimensional integral by inserting two 8- 
functions, 



The last exponential factor in this formula is nothing else but exp(i<7 ■ r) written in 
cylindrical coordinates where </> is the polar angle in momentum space while tp is the 
polar angle in coordinate space. Applying the Whittaker prescription to both sides 
of this equation we obtain on the left-hand side (up to the factor — i/kk±) the Bessel 
beam (5.8). On the right-hand side we have a superposition of plane waves whose 
wave vectors have fixed length and lie on the surface of a cone. The phases of these 
plane waves are controlled by the factor exp(iM</>). 

5.3.2. Construction of Gaussian beams from plane waves. Since Gaussian beams 
are nonmonochromatic, their expansion into plane waves will reveal their spectral 
composition. To obtain the spectrum it is sufficient to write down the function $ 
given by (5.17) as a three dimensional Fourier transform. This is done by an insertion 
of an integration over k\\ and an appropriate 8 function, 




(5.23) 




(5.24) 




(5.25) 



where we used the equality: 




(5.26) 



which is valid in the subspace defined by the S function in (5.25). This formula for ui 
determines the spectrum of the wave. Note, that for a fixed value of p the frequencies 
u> lower than c\p\ are excluded. This property is valid not only for Gaussian beams, 
but for all beams obtained from the solutions of the Schrodinger equation (5.16b). 



5.3.3. Construction of exponential beams from Bessel functions. In the previous two 
examples wc constructed solutions of the d'Alembert equation from plane waves. 
One may also construct a new solution as a superposition of other solutions, for 
example Bessel solutions. This possibility is exemplified by exponential beams [14]. 
Exponential beams are similar to Gaussian beams, but they decrease more slowly with 
the distance from the axis. They are obtained as a superposition of functions (5.7) with 
fixed values of M and k\\ but with the varying value of k±. Without loss of generality 
one may assume that A = 1 since the change of the sign of A results only in complex 
conjugation. Choosing the dimensionless integration variable as w = <Jl + (/c_l/^||) 2 
one obtains: 

/oo 
dw g(w)e- ic ^ \ wt J M (\h \ P Vw* - 1). (5.27) 

where g{w) is an arbitrary weight function. There are several choices of this function 
for which the integral can be analytically evaluated [14]. For beam-like solutions, the 
longitudinal component ku should not be spread out too much and also it must be 
much larger than the transverse component k±. The integration over w in (5.27) may 
be performed only in a few cases. The simplest choice is: 

g(w) = (w 2 - i)Af/2 e - c l fe iil™, (5.28) 

where r is a parameter with the dimension of time whose inverse determines the width 
of the beam spectrum. The formula 6.645.2 in [42] gives: 

XM „, ( ft * ,, „ = ^g^W, (5-29) 

where 

s = vV -c 2 (t-ir) 2 . (5.30) 

The MacDonald function ifjvf+i/2 f° r half-integer values of the index reduces to an 
exponential with some prefactors. In the simplest case when M = one obtains: 

Xo fc|l (p,^M)=e ifc '' z ^j^. (5.31) 

The RS vector obtained from this function by either Whittaker or Debye construction 
will exhibit an exponential fall-off with the increasing value of p. 



5.4- Construction of solutions from an arbitrary complex function 

A large class of solutions can be obtained by the Whittaker construction starting from 
the solution of the d'Alembert equation $ in the following form: 

*(r,t) = rt^-£^±^, (5.32) 




Figure 2. Lines of the electric field (left) and the magnetic field (right) 

in the y = plane drawn for the complexified Coulomb solution. The vector m 

point in the vertical direction. 



where g is an arbitrary complex function. The first term describes an outgoing 
spherical wave; the second term describes an incoming spherical wave. Both waves 
are singular at r = but their difference is regular. Perhaps the simplest solution of 
this kind is: 

$(r,t)= e - I " t . (5.33) 

r 

A more elaborate example was studied in connection with the problem of photon 
localization in [9] where the function <?(t) was chosen in the form: 



g{r) = exp(-2vTTir). (5.34) 

In Sec. 5.5.2 we give another important example of a solution defined by a very simple 
choice of g(r). 



5.5. Generation of new solutions by an imaginary shift 

Maxwell equations arc invariant under all Poincare transformations, but this obviously 
does not lead to any essentially new solutions. However, as pointed out by Trautman 
[43], all linear rclativistic equations are invariant under the complex Poincare 
transformations and this fact may be used to generate entirely new solutions. We shall 
show that even in the simplest case of an imaginary translation we obtain a different 
solution. The use of the RS vector is almost mandatory here since imaginary shift 
cannot be applied to real field vectors. We give below two examples of an imaginary 
shift. In the first case we add an imaginary space-like vector and in the second case 
it will be a time-like vector. 



5.5.1. Complexified Coulomb field. The simplest, but a highly nontrivial example of a 
new solution of Maxwell equations is obtained by an imaginary translation r r+im 
in the Coulomb field [44]. This produces the RS vector of the form: 

r + im r cos ib + m sin ib m cos tp — r sin ib 

+ i ^ , (5.35) 



(r 2 -m 2 + 2im-r) 3 / 2 R 3 R 3 



where 

R= [(r 2 -m 2 ) 2 + 4(m-r) 2 ] 1/4 , V = ^ arg (r 2 - m 2 + 2im-r) . (5.36) 

The electric and magnetic field lines as seen in figure 2 exhibit quite an intricate 
form. Since the fields possess cylindrical symmetry, the field configuration in the y = 
plane carries full information. Far from the source the electric field approaches the 
field of a point charge and the magnetic field approaches the field of a magnetic dipole. 
Since the expression (r 2 — m 2 + 2ira-r) 3 / 2 is double- valued, one must introduce a 
prescription how this expression is to be evaluated. To keep the fields continuous 
for r > \m\ one must introduce discontinuities on the disk of radius \m\ 7 as shown 
in figure 2. The discontinuity of the normal component of the electric field gives the 
surface charge 2(\m\ 2 — x 2 — y 2 )~ 3 / 2 and the discontinuity of the tangential component 
of the magnetic field gives the circular surface current 2(|m| 2 — x 2 — y 2 )~ 3 ^ 2 flowing 
in the azimuthal direction [45]. 



5.5.2. Complexified fundamental solution of the d'Alembert equation. A second 
example of a new solution generated by an imaginary shift is the one studied by 
Synge [46]. It can be obtained by the Whittaker construction from the following 
complexified solution of the d'Alembert equation: 

$(r,t)= 1 (5.37) 

&\t — lay — r z 

where a is a real constant that defines the scale of this solution. For a = this is the 
Hadamard fundamental solution of the d'Alembert equation singular on the light cone, 
but when a ^ this function is analytic in all spacetime. This function is perhaps the 
simplest solution of the d'Alembert equation with finite energy 

The electric and magnetic fields for the Synge solution are fairly complicated, but 
the energy density and the Poynting vector are given by simple formulas (c = 1): 

(a 2 - t 2 + x 2 + y 2 + z 2 ) 2 + 8t 2 (x 2 + y 2 ) + Aa 2 t 2 , 
F*-F=1Q± '- i ^ , (5.38) 



i a 



x (a 2 + t 2 + x 2 + y 2 - z 2 ) + 2azy 
y (a 2 + t 2 + x 2 + y 2 - z 2 ) - 2azx , (5.39) 
2z (x 2 + y 2 ) 

where d = ^(a 2 — t 2 + x 2 + y 2 + z 2 ) 2 + 4a 2 i 2 ^ . Note that the energy density is an 
even function of time while the Poynting vector is an odd function. The behavior 
of these functions in the z = plane explains very well the salient properties of 
the solution. In figure 3 we show the energy density for four values of time that 
characterizes a collapsing/expanding shell. This figure should be viewed from left to 
right for negative values of time ending at t = when the collapse stops (the Poynting 
vector vanishes). For positive values of time, starting with the rightmost plot, we see 
the expansion of the wave. The Poynting vector in the z = plane has no z component 
and describes the energy flowing towards the center for negative values of time and 
away from the center for positive values. The complexified fundamental solution is 
also directly related to the class of solutions described in Sec. 5.4. The solution of the 
d'Alembert equation (5.37) can be written in the form (5.32) with the function g(r) 
chosen in the form g{r) = \/{t — ia). 

It is worth noting that the solution of Maxwell equations obtained from the 
function (5.37) shows up in several unexpected places. In particular it appears in 



Figure 3. Energy density in the 2 = plane as a function of time. 



the uncertainty relation for photons [18]. Another distinct property of this solution 
is that it can be obtained from the standard plane wave by applying a conformal 
transformation [12]. 

6. Vortex lines of the electromagnetic field 

Another area of research where the description of the electromagnetic field in terms of 
complex functions is necessary is the study of vortex lines: the complex RS vector is 
indispensable in these studies. Vortex lines are intimately connected with orbital 
angular momentum and the solutions of wave equations carrying orbital angular 
momentum are most naturally described by complex functions. The signature of 
a nonvanishing orbital angular momentum is the complex factor e lMv . This factor is 
singular along a line and the removal of the singularity requires the vanishing of the 
wave function. The mere vanishing of the wave function, however, is not enough to 
create a vortex line; there must be also the circulation of velocity. 

While there is no ambiguity in nonrelativistic quantum mechanics as to how 
to define vortex lines, the problem does not have an obvious solution for multi- 
component wave functions. For the one-component wave functions in nonrelativistic 
quantum mechanics vortex lines appear generically [10]. For the multi-component 
wave functions, such as the RS vector, the requirement that all components vanish can 
be satisfied only on rare occasions. The most notable examples are Bessel beams; all 
components of the RS vector vanish for |M| > 2 along the z axis. The representation of 
Bessel beams in cylindrical coordinates (5.9) exhibits the factor e lMv> clearly indicating 
the presence of the vortex line. Other examples of vortex lines are found in null 
solutions of Maxwell equations [12]. These are the solutions with the vanishing square 
F 2 of the RS vector discussed in the next section. 

The generic vortex lines, i.e. those that are present in almost any solution 
of Maxwell equations can be obtained by replacing the quantum-mechanical wave 
function if; by the square of the RS vector [12]. Since F 2 is one complex function, 
the equations RcF 2 = and ImF 2 = define two surfaces. The intersection of these 
surfaces defines a line. The velocity field is defined by analogy with the quantum- 
mechanical case as: 

(F 2 )*d^F 2 -F 2 d, l (F 2 )* _ Sd^P-Pd^S 

W M- 2i | F 2|2 S 2 +P 2 ' 1 j 

We are able to define a genuine velocity four- vector since F 2 is a relativistic scalar field 
(2.6). As an illustration we calculate the velocity field for the RS vector given by the 
superposition of three monochromatic circularly polarized plane waves propagating 
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Figure 4. Streamlines of the velocity field for the superposition of three 
circularly polarized plane waves moving in mutually orthogonal directions. 
This pattern is periodically repeated in vertical and horizontal directions. 

ing in mutually orthogonal directions, 

F = c-' luJt [{I + im)e in k + {n + il)c im k + (m + m)e"- fc ] . (6.2) 

The three vectors I, m and n form an orthonormal set and this results in a very 
simple form of F 2 and subsequently of u M , described in detail in [11]. Vortex lines 
are aligned with the vector I + m + n. In figure 4 we show the streamlines of the 
velocity field in the plane orthogonal to this vector I + m + n. Vortex lines form a 
regular stationary lattice with the same number of vortices with clockwise and with 
anti-clockwise circulation. 



7. Knotted fields 



There is a special category of solutions of the Maxwell equations that were even given a 
special name of pure electromagnetic waves in the early days of modern electromagnetic 
theory. Nowadays they are called null waves because both invariants (2.6) are zero 
(the square of the RS vector vanishes). The simplest example of a null wave is a 
general plane wave (5.14). Bateman [27] gave a neat prescription how to construct 
a null wave in terms of the RS vector. Take a pair of complex functions a and (3 of 
spacetimc coordinates. The RS vector defined as: 

F = VaxV/3 (7.1) 

will be a null solution of the Maxwell equations if a and j3 satisfy the following 
conditions: 

Va x V/3 = i {d t aVP - d t f3Va) . (7.2) 

One may check using (7.2) that the RS vector (7.1) satisfies the Maxwell equations. 
It can be shown [48] that every null solution of the Maxwell equations can be 
represented in this form. Bateman construction is particularly useful in the description 



Figure 5. Lines of the electric field corresponding to the choice p = 2 and q = 3. 
The plot on the right shows a single trefoil artificially enhanced for clarity. In the 
plot on the left the line of the electric field was continued much further showing 
a rolling trefoil (the field line does not exactly close). 

of electromagnetic waves whose field lines possess interesting topological properties. 
Let us first note that each pair of functions a and /3 satisfying (7.2) defines a whole 
family of null solutions whose members are described by the formula 

F = h(a,p)Va x V/3, (7.3) 

where h(a,(3) is an arbitrary function. Note that we may represent every such RS 
vector again in the form (7.1) by replacing the original functions a and /3 by /(a,/3) 
and g(a, /3), where the new functions are so chosen (there are many such choices) that 
h = d a fd p g- d a gd f. 

As an example we shall consider the RS vector built according to (7.3) from the 
following functions a, (3, and h [23]: 

where p and q are natural numbers. In the simplest case, p = 1, q — 1, the field lines 
form collections of linked circles. Topologically more intricate field configurations are 
obtained when p and q are coprime. For example for p = 2 and q = 3 the field 
lines form a rolling trefoil shown in figure 5. For p — 2 and q = 5 we obtain an 
even more elaborate bundle of field lines shown in figure 6. A detailed discussion 
of these knotted field configurations can be found in [23]. Knotted electromagnetic 
field configurations not only have an interesting mathematical structure but they were 
considered as models of ball lighting [47] . 

8. The RS vector in a general coordinate system 

In this section we show that the usefulness of the RS vector extends to arbitrary 
coordinate systems (even in curved spacetime). The Maxwell equations in this case 
can also be reduced to the equations for one complex vector. 



Figure 6. Lines of the electric field corresponding to the choice p = 2 and q = 5. 
The plot shows the bundle of electric field viewed from the top (along the z-axis) . 



The original Maxwell equations (1.2a-d) hold in any coordinate system also in 
curved spacetime. The nabla operator that appears in these equations denotes the 
derivatives with respect to the chosen coordinates and the field vector have components 
along the coordinates. For example, in spherical coordinates the nabla stands for 
(d r , dg,d v ). The whole information about the coordinate system is in the constitutive 
relations between the (D, B) pair and the (E, H) pair. The constitutive relations (cf. 
for example [5, 8, 49]) involve the components of the metric tensor g^ v and its inverse 
g^ v (the Greek indices run from to 3), 

1 ( , 



Ei = { 9ij \~— + '-fro,,!! 1 ) , (8.1a) 
H < = 4o (-9ii/V=9— - '■<i ,l <,, l l) J ) , (8.16) 



.a 00 V " - n 

where wc used the summation convention over repeated indices (Latin indices run from 
1 to 3) and g is the determinant of the g^ v matrix. The space with a nonvanishing 
space-time component <?ofc of the metric tensor is bi-anisotropic because the electric and 
magnetic components are mutually coupled through the constitutive relations. The 
vector components in these formulas are the coordinate components not the physical 
components defined as projections on an orthonormal set of basis vectors that were 
used so far. 

The expressions (8.1a-b) can also be rewritten in terms of complex vectors J- 
and Q built according to the formulas (1.1) and (1.4), but with the Cartesian vectors 
(D, B) and (E, H) replaced by their counterparts in a general coordinate system, 

Qi = SijP = {-ga/y/^g + ig ok e kij ) T\ (8.2a) 

r = (S-'fQj = — {-g lJ V~9 - W Hj ) Q r (8.26) 

We use new symbols to denote these vectors because they should not be identified with 
their Cartesian counterparts. Using Cartesian coordinates we have not distinguished 



between covariant and contravariant components; this distinction is necessary in 
curvilinear coordinates. 

The easiest way to find the relations (8.1a-b) is to use the variational principle. 
The Lagrangian density in spacetimc endowed with the metric tensor g^ v has the form 

c = -iV^ag^^Uufxp, (8.3) 

where f^ v is built from the four- vector potential / M „ = d^A v — d„A^. The vectors 
(E,B) are obtained from the components of / M „ as follows: E k = fok, B k = —fy. 
The form of the Lagrangian (8.3) is a unique generalization of the expression in the 
Cartesian coordinates to an arbitrary coordinate system. The Maxwell equations in 
the tensor form derived from this Lagrangian are: 

d^fux + dxU„ + d x f^ = 0, (8.4a) 

dphr = 0, (8.46) 
8C 

hflv = ~ 2 ^r- = V=$ g^g up h P • (8.4c) 



From these equations we can read of the relations (8.1a-b), after the identification: 
D k = h k0 and H k = -h ij . 

The matrix that determines the constitutive relations (8.2a-b) appears also 
in Maxwell equations for the RS vector. These equations may be expressed either in 
terms of J- or in terms of £?, 

id t r = e ijk dj£ k iP, (8.5a) 

\at(E-*fGj = e ijk djg k . (8.56) 
In each case one has only one complex vector to describe the electromagnetic field as 
was announced at the beginning of this section. 

The definition of the RS vector in a general coordinate system involves a subtle 
point that should be elucidated. We shall start by writing the formula for the 
Hamiltonian in a general coordinate system, 

H = \ J d 3 r(D-E + B H) = \ J dx 1 dx 2 dx 3 (D i E i + B i H i ) 



= / dx l dx 2 dx 3 F n £ l3 F J , (8.6) 



where we used the constitutive relations (8.1a-b). The Hamiltonian coincides with the 
energy when 500 = 1 and all components f/ofc vanish. These conditions are satisfied 
in the simple case of fiat space that will be considered first. The subtle point is 
the presence of different volume elements in the expressions (2.5a) and (8.6). In the 
first case it is d 3 r; in the second case it is dx 1 dx 2 dx 3 . After the transformation to 
the coordinates x l the volume element d 3 r contains in addition to dx 1 dx 2 dx 3 also 
the Jacobian of the transformation. For example, in spherical coordinates we have 
d 3 r = r 2 sin# drdddip. The Jacobian does not appear explicitly in the expression (8.6) 
for the energy, it is contained in the product of the RS vectors. 

8.1. Bessel beams in cylindrical coordinates 

We shall apply now this formulation to the cylindrical coordinate system in flat space. 
The time variable in this case is not affected and the spatial line element is: 

ds 2 = dp 2 + p 2 dip 2 + dz 2 . (8.7) 



Hence, \J—g = p and the constitutive relations read: 

( T p T z 
{G P ,G v ,Gz) = — ,p^, 



P P 



The equations for T become: 



i-d t T p = d v — - d,pFf, (8.9a) 
c p 

1 T p T z 

i-d t f = d z 8 P — , (8.9b) 

c p p 

i-d t T z = dppF? - d v — , . (8.9c) 
c p 

Our aim now is to reproduce the results for the Bessel beam. To this end we shall 
seek the solution of equations (8.9a-c) in the form: 

F(p t <p,z,t) =e iX ^ z - u,t+M,fi) T{p). (8.10) 

Upon substituting this Ansatz into (8.9a-c), we obtain: 

\M 

kT p = — J" 2 - ihpT v , (8.11a) 
p 

\tpt> = llp» _ \d p -T z , (8.116) 
P P 

kT z = XdopJ* - — T p . (8.11c) 
P 

Next, we solve the first two equations (8.11a) and (8.116) with respect to T p and T v , 

TP = TT ( kM + X P k \\ d p) — > (8.12a) 
kj_ p 

pF v = -^T {k\\M + Xkpd p ) —. (8.126) 
fe_i_ p 

After the substitution of these expression into the third equation we obtain the Bessel 
equation for the function T z / p, 

1 M 2 \ T z 

-d pP d p - -=- + k\ — = 0. (8.13) 
P P J P 

Comparing this with the solution for the Bessel beam (5.9) one can see that T z in 
cylindrical coordinates differs by a factor of p from the Cartesian F z . By comparing 
the expression for the energy (2.5a) and (8.6) one can establish the full connection 
between T and F, 

[F p , J", T z ) = (pF p , F v , P F Z ) . (8.14) 

8.2. Multipole fields in spherical coordinates 

The spatial line clement in spherical coordinates is: 

ds 2 = dr 2 +r 2 dif 2 + r 2 sm6 2 dip 2 . (8.15) 
Hence, y/—g = r 2 sin and the constitutive relations read: 

(G r ,Ge,G,p) = f-5?-T,^:,sin^V (8.16) 



The equations for J- after the separation of the variables t and if became: 



kT r = dgsiaej^ - 



sin 9 



kP = ^r^ -drwnej^, 



kT v = a 



-do 



(8.17a) 
(8.176) 
(8.17c) 



sin a r^ sin t) 

These equations are satisfied by the expressions (5.12) provided one takes into account 
the following relations between T and F: 

(F r ,F d ,F v ) = (r 2 sin 6F r ,r sin eF ,rF v ) , (8.18) 

which may be obtained from the two formulas for the energy as in the derivation of 
(8.14). 



8.3. Electromagnetic field in curved space 

We consider here, as a simple example, the curved space in the vicinity of a black hole 
of mass M. We use the standard form of the Schwarzschild metric, 



2 i 2 

c dr 



c 2 (l-r s /r)dt 2 



dr 2 



(d9 2 



sin 2 9 dip 2 ) 



(8.19) 



1 - r s /r 

where r s = 2GM/c 2 is the Schwarzschild radius. The Maxwell equations in this case 
are: 



d t T r = dg sin 9 (l-— ) 



1 - — 

r 



sint 



-d t F u = a 



v r 2 sin 9 



-a.sin6» 1- — )T lp , 



aj" = a. i - — 



r 

r 



-Wujt 



Y lm (6,ip). 



(8.21) 



(8.20a) 
(8.206) 

de^— n - (8.20c) 
sm U r z sm 9 

We shall seek the solutions of these equations in the spherical form (patterned after 
equations (5.12) and (8.18)): 

-Z(Z + l)rsin0/ r (r) 
-Xkrdtpf^r) - sin 6* d g d r rfe(r) 
Xkrdef^r) - 1/ sva.6 d v d r r fe(r) 

After cumbersome calculations we obtain three equations for the radial functions 

{fr, fe, ftp), 

rfr(r)-(r-r s )U(r)=Q, (8.22a) 
[rfe(r)}'-[(r-r s )U(r)]' = 0, (8.226) 
r 2 (r - r s )f' e f (r) + (2r 2 - rr s )f' e {r) + r s f e (r) - 1(1 + l)r/ r (r) + fcV/^r) = 0, (8.22c) 

where primes denote derivatives with respect to r. After the elimination of fe and f v , 
we obtain the following equation for the function g = rf r : 

k 2 r 2 1(1 + 1) 
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This equation, known as the confluent Heun equation, does not have solutions in terms 
of known functions. 
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Figure 7. Effective potential V (measured in units of rj 2 ) as a 
function of the distance from the black hole (measured in units of r s ). 
The potential curves are drawn for three values of I. 



The physical meaning of the solutions may be best seen after the change of 
variables r* = r + r s ln(r/r s — 1). In the new radial coordinate r* the horizon, that lies 
at r = r s in the standard Schwarzschild coordinate, is pushed to — oo. The function 
r*(r) cannot be explicitly inverted but r(r*) can be found numerically 

The equation in the variable r* has the form of the one-dimensional Schrddinger 
equation [50], 



where r is to be treated as a function of r* and g(r*) = g(r(r*)). The second term on 
the left-hand side plays the role of the potential V(r*). The function V(r*) plotted 
in figure 7 has the form of a potential barrier. The maximum of the potential barrier 
has the value 41(1 + l)/27 and is located at r* = (3/2 - ln2)r s w 0.8r s . Due to the 
repulsive character of the potential some part of the electromagnetic waves impinging 
on the black hole will bounce off; the rest will be sucked into the black hole. The 
reflection by the potential barrier is more pronounced for small values of k and for 
large values of I. The insights drawn from the analogy with the Schrddinger equation 
are fully justified because the RS vector (represented here by g(r*)) takes on complex 
values so that it can be treated as an analog of the Schrddinger wave function. 

9. The RS vector in a relativistic setting 

The three-dimensional RS vector is a part of the relativistic antisymmetric complex 
tensor F^" , 



where e A " yAp is the Levi-Civita symbol (antisymmetric in all indices and e 0123 = l). 
The Levi-Civita symbol by itself is not a tensor but a tensor density [50]; it becomes 
a regular tensor after the division by \J—g. Therefore, to be precise, F yiv is the tensor 
density; it is a sum of two tensors multiplied by \J—g. The definition of is valid 
in an arbitrary coordinate system, also in curved space. 




(8.24) 




(9.1) 



Since the f\ p has 6 real components, only three components of the complex tensor 
F^ v can be independent. Indeed, this tensor is self-dual; it transforms into itself under 
the following duality transformation: 

^\F xp = F^. (9.2) 



2x^5 Ap 

The relativistic tensor F^ v obeys Maxwell equations in the form: 

dpFi™ = 0. (9.3) 

The real and the imaginary part of this complex equation reproduce the two pairs of 
Maxwell equations written in the relativistic notation, 

^ = 0, (9.4) 

d„fux + dxU„ + = °- ( 9 -5) 

The three components of the RS vector can be identified with the space-time 
components of F liV , 

F=(F W 7 F 20 ,F 3Q ). (9.6) 

However, due to the condition of self-duality (9.2), 

/p23 p3l pl2\ = i tplO F 20 5 F 30) /g ?) 

the components of the RS vector can also be identified with the space-space 
components of F^ . This explains why the components of the RS vector under Lorentz 
transformation transform among themselves (2.8), even though they are at the same 
time the components of a relativistic tensor F^ IU whose space-time components should 
mix with space-space components. In other words, the electric and the magnetic field 
vectors can be combined into two equivalent geometric objects: a complex vector F 
transforming under the complex orthogonal group 0(3, C) or a real second rank tensor 
fp U transforming under the real group 0(1, 3). 

In flat space and in Cartesian coordinates the Fourier representation (3.7) written 
in term of F^" has the following explicitly relativistic form [6]: 

F^{r,t) = jA^J ^feT{k) [/ + (fc)e ife — i " rf + /* (fc)e-** +1 * rt ] . (9.8) 

The tensor e^(k) is self-dual; its components can be expressed in terms of the 
polarization vector e(k). 

Quantum theory 



10. Canonical quantization of the electromagnetic field 

The standard procedure to formulate the quantum theory is the canonical quantization 
of the classical theory. This procedure was first applied to electromagnctism by 
Hciscnbcrg and Pauli [51]. Canonical quantization requires the identification of the 
canonical variables and this was properly done in the general case by Born and 
Infold [52]. The Maxwell equations clearly indicate the right identification: canonical 
variables are the field vectors D and B whose time derivatives are determined by the 
evolution equations. It does not really matter whether one assigns the role of position 



to B and the role of momentum to D or vice versa. However, the first choice is more 
natural since the sources — the counterparts of the external mechanical forces — appear 
in the equation for D. Either way, what is only needed are the Poisson brackets for 
canonical variables that will enable one to determine the commutators of the field 
operators. The basic requirement for the Poisson brackets is that they reproduce the 
equations of motion, 

d t Q(t)={Q(t),H}, (10.1) 

where Q(t) is any variable and H is the Hamiltonian of the system. The Poisson 
bracket for B and D in the form: 

{B i (r,t),D j (r , ,t)} = e ijk d k 5^(r - r') (10.2) 

not only reproduces the time dependence (10.1), but it also assures the proper action of 
the remaining generators of the Poincare group. Note that the Poisson brackets (10.2) 
are universal — they contain no material constants. This validates the choice of B and 
D (and not, for example, H and E) as the canonical variables. According to the rules 
of canonical quantization, the Poisson brackets are replaced by the commutators of 
the corresponding operators (divided by ih). The commutation relations for B and 
D lead to the following commutators for the operators F and representing the RS 
vector and its Hermitian conjugate: 

[Fi(r,t),F](r',t)} = \tvce ijk d k 5{r - r'). (10.3) 

The photon picture emerges more directly when one uses the RS vector than in the 
standard formulation. After the Fourier transformation of the RS operators according 
to the formula (3.7) we arrive at the expression, 

F(r,t) = J J0j^e(k) [f + (k)c ik - r - +/t(fe) e -*-+-*] , (10.4) 

where the classical amplitudes arc replaced by the operators. In order to recover 
the commutation relations (10.3) for the RS operator one must impose the following 
commutation relations on the operators fx- 

[f\(k),fl(k')} = faj5 xxl S^(k-k'). (10.5) 

Were it not for the factor tuo these would be exactly the commonly used commutation 
relations for the photon annihilation and creation operators. We shall get rid of the 
factor fie by the following rescaling: 

A(fc) - Vh~ca x (k), f*x(k) = yfhca{{k). (10.6) 

The commutation relations for the annihilation and creation operators obtained after 
the rescaling are: 

[a x {k),a\,{k')\ = k6 xy S^(k - k'), (10.7) 

We retained the factor of k in these relations because its presence is connected 
with the relativistic normalization of the photon wave functions, as will be shown 
in Sec. 11.1. The formulas (10.6) justify our notation in (3.6) where we introduced 
complex conjugation in the second term of the Fourier transform. Finally, the Fourier 
representation of the RS field operator takes on the following form: 

, r r! 3 k r i n 

F(r,t) = Vtw — ^e(fe) a + (fc)e ifc ^ iwt + a T _(fc)e- ifcr+i " t . (10.8) 
J (27r) cS /^ L J 



The time dependence of the annihilation and creation operators is given by exp(— iut) 
and exp(iwi), as it should be. Quite often this time dependence is used to distinguish 
the annihilation from the creation operators. However, this method sometimes fails 
while the assignment based on the canonical commutation relations always works [16]. 

The generators of the Poincare transformations in quantum electrodynamics are 
obtained by replacing the classical RS vectors in (2.5a-d) by quantum field operators. 
Expressed in terms of the creation and annihilation operators, the generators have the 
form: 

& = E I —a{(k)hujax{k), (10.9a) 
a[{k)hka x (k), (10.96) 
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fk 
d 3 k 
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M = ^ / — a\(k) [iHD\ xk + Xhn] a\(k), (10.9c) 

A •* 

Cr1 3 k 

™ = J2 -ra\{k)\hwD x a x {k). (10.9d) 

A •* 

The normal ordering of the creation and annihilation operators removes the unphysical 
vacuum contributions. 



11. Quantum mechanics of photons 

Quantum theory of the electromagnetic field was obtained in the previous section 
directly from the classical Maxwell theory and not by the second quantization of 
quantum mechanics of photons. Usually, second quantization upgrades the wave 
function to the status of the field operator. However, the photon wave function is 
a controversial notion because it is impossible to endow it with all properties of wave 
functions for massive particles. Nevertheless it is possible to derive unambiguously 
quantum mechanics of photons by "reverse engineering" from quantum field theory. 
The RS vector is very useful in this reconstruction. To proceed, we consider the one- 
photon sector of Hilbcrt space and we will define the action of the important operators 
in this sector. 



11.1. Photon wave function in momentum representation 

The most general one-photon state |/) is generated from the vacuum state |0) by the 
action of a superposition at of creation operators, 

— / A (fcK(fc)io). (ii.i) 

The normalization of the state vector |/) translates into the following normalization 
condition for the functions f\(k): 
_ r r fih 

ii/ii 2 = E ^imh(k) = i. (n.2) 

A J 

We used the same symbols f\(k) as in the Fourier expansion of the RS vector (3.7) 
because the two-component function f\(k) determines the state of the photon and 
fully deserves the name of the photon wave functions in momentum representation. 



The existence of the photon wave function in momentum representation is a direct 
consequence of the theory of representations of the Poincare group and its subgroup: 
the Lorentz group. According to this theory, fully developed by Wigner [53], there 
are two one-dimensional representations f±(k) corresponding to two circular photon 
polarizations. Under a Poincare transformation IT the transformed functions ' f± are 
multiplied by fc-dcpendent phase factors, 

7±(nfe) - cx P [±i0(n,fc)]/±(fc). (ii.3) 

These functions should be identified with the Fourier amplitudes of the RS vector 
appearing in (3.7) because they have identical transformation properties under the 
Poincare group. However, we shall not prove it here since to prove that the 
transformation properties of the electromagnetic field under the Poincare group induce 
appropriate phase transformations of the Fourier amplitudes is a rather cumbersome 
task. The transformation properties of the annihilation and creation operators must 
be the same as those of the corresponding functions / since the RS operator must 
transform in the same way as its classical counterpart. 

The expectation values ((f\H\f), etc.) of the quantum generators (10.9a-d) in 
one-photon states coincide with the classical expressions (3.9a-d). At this point we 
would like to comment on the significance of the factor 1/k in the volume element in 

(11.1) . This factor produces an invariant volume element on the light cone and the 
whole formulation is relativistically covariant. The expression f^(k)f\(k) d 3 k/k is the 
invariant probability to find the photon with momentum hk and polarization A in the 
volume d 3 k provided the function f\(k) is normalized. The normalization condition 

(11.2) is relativistically invariant, as it must be. 

The action of the generators of the Poincare transformations (10.9a-d) in 
quantum electrodynamics on the state vector |/) is equivalent to the following action 



on the functions f\{k): 

Hf x (k) =tiufx(k), (11.4a) 

Pf x (k) =hkf x (k), (11.46) 

Mfx(k) =H\kx jDx + Xn^j / A (fc), (11.4c) 

Nfx(k) =hwiD x fx{k). (11.4d) 



We would have obtained the same results from the Wigner theory of representations. 
This makes the whole scheme fully consistent and one may identify the amplitudes of 
one-photon states (11.1) with the photon wave functions in momentum space. Note 
that the quantum theory of the electromagnetic field could be, in principle, obtained 
by second quantization of the quantum mechanics of photons, i.e. by upgrading 
the photon wave functions in momentum space f± to the status of the annihilation 
operators. 

11.2. Photon wave junction in coordinate representation 

The definition of the photon wave function in coordinate representation has well 
known difficulties. The straightforward Fourier transformation of the momentum 
wave function does not produce a field that transforms locally under the action of 
the Poincare group. In order to identify the mathematical function that may serve as 
the photon wave function in the coordinate representation, one has to decide which 
properties one wants to keep and which properties one is willing to abandon. In 



our opinion a natural candidate is the RS vector. This decision is supported by the 
analogies between the description of the electromagnetic field in terms of the RS vector 
and quantum mechanics. 

Of course, Silberstein writing his paper twenty years before the appearance of 
the wave function, could not have noticed the mathematical analogy with quantum 
mechanics, although his vector has many features of the quantum-mechanical wave 
function. First of all, the RS vector satisfies the evolution equation (1.6) that has the 
form of the Schrodinger equation with the Hamiltonian H = p-S, 



ihd t F(r,t) = - V-SF{r,t), 

1 

where S are the spin-one counterparts of the Pauli matrices, 
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The analogy between the scalar product in quantum mechanics and the expression 

(F 1 \F 2 )= Jd 3 rF*(r)-F 2 (r) (11.7) 

enables one to introduce such notions as the complete set of solutions and change of 
basis. The most important change of basis is the Fourier transformation. When such a 
transformation is applied separately to the real electric and magnetic held vectors, the 
Fourier components are not independent and one has to impose cumbersome reality 
conditions. This problem does not exist for the RS vector. The only constraint that 
has to be imposed is the transversality of the Fourier transform: k-F(k) = 0. 

The analogy between the probability density ip*ip and F*-F can also be illustrated 
with the phenomena of wave packet spreading during the free time evolution. It 
follows from the standard Schrodinger equation that the dispersion of the probability 
distribution a 2 = ((r — (r)) 2 } has the following time dependence [54]: 

o 2 r = a 2 \ t=tl + (t - hful/m 2 , (11.8) 

where a 2 is the dispersion in momentum which for the free evolution is a constant 
of motion and t\ is the time at which the wave packet attains its minimal size. 
An analogous formula holds in electromagnetism since the second derivative of the 
dispersion is proportional to the total energy. It follows from the Maxwell equation 
(1.6) that 

^ J d 3 r{r - (r(t))) 2 F*(r,t) ■ F(r,t) = c 2 J ' d 3 r F* (r,t) ■ F(r,t). 

This formula shows a quadratic time dependence because the total energy is a constant 
of motion. Thus, as in quantum mechanics, one obtains: 



(11.9) 



J d 3 r(r- (r(t))) 2 F*(r,t) ■ F(r,t) 



d^ir-irih^F^r,^) 



F(r, h) + (t- hfc 2 J d 3 r F*(r, 1) ■ F(r, t), (11.10) 

where ti is the time at which the dispersion attains its lowest value. 

In spite of these analogies there are important differences between the Schrodinger 
wave function ip and the RS vector. Even at first glance one can see that they differ in 
having different dimensions: F*(r, t)-F(r, t) is the energy density while V'*( r i t)4 ! { r i t) 
is the probability density to find the particle at the position r. 



A more important difference stems from the nonexistence of a bona fide photon 
position operator. In particular, a natural choice — an operator r whose action 
on F(r,t) results in the multiplication by r as in quantum mechanics — is ruled 
out since the multiplication by r produces a function that is not divergenceless: 
V • XiF(r,t) = Fi(r,t). The nonexistence of the position operator precludes the 
existence of the photon probability density. Nevertheless, one may formulate the 
uncertainty relation for photons using the energy density. 

11.3. Uncertainty relation for photons 

In most general terms, the uncertainty relation for position and momentum imposes 
limitations on the spread of quantum systems both in the coordinate space and in 
the momentum space. In its standard form, this uncertainty relation involves the 
dispersions in position and momentum, but other measures of the spread (for example, 
cntropic measures [21]) have also been successfully used. Wc proposed in [18] to 
base the uncertainty in position for photons on the energy density — the modulus 
squared of the RS vector. More precisely, we introduced the second moment of the 
energy distribution as a measure of the spread of the photon extension in coordinate 
space. The corresponding quantity is denoted by Ar (and not Ar 2 ) because it has the 
dimension of length, 

Ar= ^lW / d3rr2 - F *( r ' < )- F ( r ' i )- ( 1L11 ) 

Division by the norm is necessary to make Ar independent of the overall amplitude of 
the RS vector. In other words, Ar is normalized per one photon. A natural partner 
of Ar in the uncertainty relation is Ap, 

i r ^3 ;„ 

A P= j^pE J—h\k-(k)\\f x (k)\ 2 7 (11.12) 

that measures the spread of the photon wave function in momentum space. The lower 
bound for ArAp was evaluated analytically in [18] in the center of momentum frame 
((fc) = 0) and it gives the uncertainty relation in the form: 

ArAp > Ati. (11.13) 

The lower bound is much larger then 3/2h appearing in the three-dimensional 
Heisenberg uncertainty relation in nonrelativistic quantum mechanics. This is due 
to the fact that photon owing to its relativistic nature hinders localization. The 
photon state that saturates this inequality is the Synge solution of Maxwell equations 
[18] described in Sec. 5. 

12. The RS vector and the spinor calculus 

In order to underscore the interpretation of the RS vector as the photon wave function 
we shall now prove that it fits very well into the general scheme of wave functions 
for masslcss relativistic particles. Relativistic properties of fields and particles find 
their apt description in terms of spinors and the RS vector does not depart from 
this general rule. The RS vector is directly related to the second rank relativistic 
symmetric spinor field (j}AB{f,t). Such a spinor field has only three independent 



components (0oo, 0n, 0oi)- Assuming the relations between 4>ab and the RS vector 
in the following form: 

F x = 0n - </>oo, ^ = -1(011 + 000), ^ = 2001 = 2010, (12.1) 

we can prove the equivalence of the spinorial equations: 

d t (t>AB = -c(<T-V)f<t>cB, (12.2) 

with the Maxwell equations. To prove this equivalence we rewrite (12.2) in full matrix 
form: 



000 


001 






= — c 


001 


011 





,(12.3) 



<9 2 0oo + ( d x ~ i<9. y )0oi <9 2 O i + (d x - i<9 y )0n 
(d x + i<9 y )0oo - 5 2 0oi (d x + i<9,,)0oi - 9 z 0n 

where we employed the standard representation of Pauli matrices. Next, we use the 
relations (12.1) to reproduce the Maxwell equations (1.6). The symmetry of the matrix 
on the left-hand side of (12.3) imposes the condition that the right-hand side is also 
symmetric, i.e, 

9 2 0oi + (d x ~ id y )tj>u = (d x + i<9«)0oo - <9 2 0oi- (12.4) 

This equality implies that the divergence of the RS vector vanishes. Thus, there is 
a full equivalence of the description of the electromagnetic field in terms of the RS 
vector and in terms of the symmetric second rank spinor. However, when it comes 
to concrete applications the use of the RS vector is much more convenient due to its 
direct connection with the electric and magnetic field vectors. 

There is a simple connection between the spinor <\>ab and the complex tensor F^ v 
defined in (9.1) that can be expressed in terms of Pauli matrices combined into the 
spin-tensor S^ vAB , 

S ak = -S k0 = -ia 2 a k , SP = -St* = «' ; W A . (12.5) 

where we suppressed the AB indices. With the use of (12.1) we obtain: 

F» v = S^ AB <j>AB- (12.6) 

Null electromagnetic waves have a very simple form in terms of spinors. Namely, they 
are described by just one first rank spinor field 4>a (0ab = 0A0s)- Of course, not 
every spinor field 0a gives rise to a null solution of Maxwell equations. The second 
rank spinor 4>ab must still satisfy (12.2). 

The equation (12.2) is a special case of the general wave equation for a massless 
particle. Regardless of the spin value s, these equations have always the same general 
form. 

dt<j>AB...w = -c(a ■ V) A Z 4>zb...w- (12-7) 

They differ only in the rank of the spinor that is equal to 2s. The first-rank spinor 
describing a spin 1/2 massless neutrino satisfies the Weyl equation 

d t <t> A = -c(a- V) A Z <t> z . (12.8) 
The fourth-rank spinor satisfies the equation: 

d t <t>ABCD = -c(a ■ V) A Z 4>zbcd- (12-9) 
It describes gravitons — the quanta of the gravitational field. 



13. Summary 



In this review we have shown that the RS vector is an ideal tool to unify the classical 
Maxwell theory and the quantum mechanics of photons. On the one hand the RS 
vector gives a concise description of the classical electromagnetic field. On the other 
hand it is the best possible photon wave function [5]. In the study of solutions of 
the Maxwell equations complex functions to describe real fields are widely used since 
"Exponentials are much easier to manipulate than sines and cosines" [55]. However, 
in the standard approach one only keeps the real part as physically significant. The 
RS vector offers all advantages of working with complex functions but there is no 
redundance in this description; the real part and the imaginary part are both physical. 
Nevertheless, if the role of the RS vector were restricted only to classical theory, we 
would have to agree with Stratton [38] who in his monograph wrote: "It has been 
shown by Silberstcin, Batcman, and others that the equations satisfied by the fields 
and potentials may be reduced to a particularly compact form by the construction of 
a complex vector whose real and imaginary parts are formed from the vectors defining 
the magnetic and electric fields. The procedure has no apparent physical significance 
but frequently facilitates analysis" . In our review we have argued that it is the role 
of the RS vector as the photon wave function that brings forth physical significance 
denied by Stratton who was not aware of a dual role of the RS vector. The use of 
the RS vector as a photon wave function enables one to incorporate the quantum 
mechanics of photons into the general scheme valid for all relativistic particles. 
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